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1. SUMMARY 

During this reporting period, slgniflcunt progress has been made 
in our understanding of the web growth process. Thermal models have 
been developed that accurately predict the thermally generated stiessea 
in the web crystal which, if too high, cause the crystal to degenerate. 
The application of the modeling results to the design of low-stress 
experimental growth configurations will allow the growth of wider web 
crystals at higher growth velocities. 

A new experimental web growth machine was constructed, this 
facility includes all the features necessary for carrying out growth 
experiments under steady state thermal conditions, Programmed growth 
initiation has been developed to give reproducible crystal starts. 

Width control permits the growth of long ribbons at constant width. 

Melt level is controlled to 0.1 mm or better. Thus, the capability 
exists to grow long vreb crystals of constant width and thicltness with 
little operator intervention, and web growth experiments can now bo 
performed with growth variables controlled to a degree not previously 
possible. 
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2. INmODUCTION 


Silicon dondrlclc web la a single crysdal silicon ribbon 
mauerial which provides substanuial advantages for low-cost manufacture 
of solar cells. A significant feature of the process is the growth from 
a melt of silicon without constraining dies, resulting in an oriented 
single crystal ribbon having excellent surface features. In common with 
other more classical processes such as Czochralskl growth, impurity 
rejection iirto the melt permits the use of less pure "solar grade" 
starting material without significantly affecting cell performance. A 
unique property of the dendritic web process is the growth of long 
ribbons of controllable width and thickness which not only facilitates 
automation of subsequent processing into solar cells, but also results 
in high material utilization since cutting and polishing are not 
required. 

During the previous program (DOE/JDL Contract No. 954654), most 
of the component elements for the reproducible and steady state growth 
of high-quality web crystals were developed and demonstrated. Area 
throughputs greater than 25 cm^/mln were demonstrated for short periods 
of time. Molt replenishment for periods of up to 17 hours (a one-day 
growth cycle) was demonstrated. Thermal models were developed for 
calc\»lating temperature distributions in the web crystal as a function 
of configuration parameters. 

On the present contract and during this reporting period, three 
broad areas uf work are emphasized; 

1. The development of thermal stress models in order to understand 
the detailed parameters which generate buckling stresses. The 
model can th^jii;, be used to guide the design of Improved low- 
stress web growth configurations for experimental testing. 
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2. Experiments to Increose our understanding of the effects of 
various parameters on the web growth process. 

3. The construction of on experimental web growth inaehine which 
contains in a single unit all the mechanical and electronic 
features developed previously so that experiments can be carried 
out under tightly controlled conditions. 

Thus, the principal objective of this work has been to expand 
our knowledge and understanding of both the theoretical and experimental 
aspects of the web growth process to provide a solid base for substan- 
tial Improvements In both area throughput and web crystal quality. 
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3. 


PROGRESS IN WEB GROWTH RESEARCH 


3*1 Develqpraent and Applieation of Thetraal Stress Mod eling 
3.1*1 Basic Theory 

Although surface tension forces may be the ultimate factor that 
limits the width of dendritic web crystals, thermally generated 
stresses are presently a much more serious limitation. As the web 
crystal widens, the magnitude of the stress Increases while at the same 
time the "stiffness" of the ribbon decreases until at some critlcol 
point the growing crystal buckles. Once the ClTySbiaX liw0 deformed, it is 
not only unsuitable for use in device fabrication, but also will lose 
its single crystal character. The overall goal of the thermal stress 
modeling work in this program is to reduce the thermal stresses in the 
growing web by changing the temperature profile in the ribbon through 
use of an appropriate thermal environment. 

Thermal stressco have two major effects on the growth of ribbon 
crystals such as web: residual stresses whit*h remain in the grown 
crystal and the deformation or "buckling" which occurs during growth. 

In practice, the residual stresses can be controlled by tailoring the 
thermal environment near the growth front by proper design of the slot 
configuration in che lid and through control of the melt height. 

Even when residual stresses are not a problem, buckling still 
occurs witen the web has reached some critical combination of width and 
thickness. So far, control of the stresses causing buckling has not 
been obvious. Our goal in the present program is to identify the 
critical ouckling stress magnitude, the temperature profile necessary to 
avoid deformation, and finally the growth configuration required to 
generate this profile. 


4 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Unlike the tCBldual attesa resulting from plaatic deformation of 
the ribbon» the atrooaeo causing buckling are affected by the tempera** 
ture diotribution away from the interface region* R. W. Gurtler^^^ hoe 
modeled the atrcaaos from two temperature dlotrlbutlono which were 
identical near the melt surface but different at greater dlstoncca* The 
first distribution decreased linearly along the length of the ribbon; 
there woo no elastic stress and there would be no buckling. Unfortu- 
nately, this temperature distribution is unrealistic since absolute zero 
is reached a few centimeters from the growth front. The temperature 
decline could not bo mode more gradual without eKcesolvcly slowing the 
growth rate. 

The second temperature distribution was identical to the first 
near the growth front but parabollcally approached ambient temperature 
instead of linearly reaching absolute zero; this temperature profile 
caused large clastic stresses. 

Of course, the mere existence of thermal stresses In a growing 
web crystal does not automatically mean that the web will have residual 
stresses or deform during growth. In order for residual stress to be 
observed, the thermal stress must have exceeded the yield stress; for 
deformation to occur, the thermal p^'rosses must exceed some critical 
buckling stress. ^ determination of the buckling stress criterion and 
the required temperature distribution would be a valuable guide in the 
design of high throughput growth configurations. The following sections 
review the mathematical formulation of both thermal stress generation 
and the ribbon deformation, and elucidate several approaches to solving 
the problems. 

Stress model . The stresses in a thin ribbon of thickness t, 
length coordinate x, and width coordinate y, are related to the plane 
temperature distribution T (x,y) by means of the Airy stress function <|); 
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r 4 

0 M 

( 1 ) 

V ■ 

whece (ijjjj is the longitudinal atresfl, 0^y is the aheat atreaa, ond Oyy 
the tcanaverse a tress. These stresses are doterrained 0rom the tempera- 
ture diatributlon via the differential equation 

*. «.d E AT (2) 


where « is the thermal expansion coefficient, e is the modulus of 
elasticity (Young's modulus), 

^2 a2 

A ~ 3 a. 3 a. 3 

(I s “r T 

3x^ 9/"^ Bz^ 

is the Laplaclan or horroonic operator and A^ is its square or the 
blharraonlc operator. Boundary conditions require that both <|) and its 
normal derivative should vanish on the growth front and on the edges of 
the web. 

With proper shield design the web temperature should be nearly 
constant in the crosswise, or y, direction. For this case, T ° T(x), 
Boley and Miener^^^ have found an approximate series solution for 
equation 2 which is valid away from the melt surface. Numerical 
solutions can also be obtained by various techniques; for example, the 
Westlnghousc WECAN finite element code has been used to solve for the 
stress distribution both near at,d far from the growth front with an 
arbitrary temperature distribution for input. 

Buckling Model . Whether or not the stresses determined by 
equation 2 would cause buckling depends on another differential equation 
which also involves the biharmonic operator^^^ ; 
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whore V l8 ?oio8on*o ratio, E Id Young^s modulus, and w(x,y) is the trans” 
voiae dofleetlon (z dirocCion) of a buckled thin plate of thlcknodd t. 

The syatem of aquations 0, 7, vand 8 under the appropriate boundary 
GondltionR form the raatherautical model of buckling in silicon web growth. 


Since the web edges consist of dendritos which are loss flexible 
than the web itself, the closest model for the boundary conditions would 
seem to be edges supported by elastic bearasC^). 



± D 


.3 J 

ti W * tJ w 

Dy Ox Oy 





(4) 


on the edges y ® » c, where El is the flexural rigidity of the dendrite, 
A is the cross-sectional area of the dendrite, and D » ET^/12(1“V^) is 
the flexural rigidity of the plate or web. Here wo have added the last 
terra on the right side of equation 4 to account for the longitudinal 
variation of o^jx* Equation 4 accounts for the twisting moments and the 
vertical shearing forces at the edges. The bonding raoraonus are 
considered in an additional boundary condition; 



-s2 


2 2 
» D + V — y) 
8x^ By 


(5) 


on y “ A G, where C is the torsional rigidity of the dendrite. 

For the edge of the initial button, boundary conditions similar 
to equations 4 and 5 would be appropriate when the growth has recently 
begun. After the web is long enough to wrap around the pickup reel, the 
leading edge of the web would be fixed and the boundary conditions would 
then be 

w » 3w/Bx “ 0 
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Finally, at the growth front x - 0, the edge Is free and the boundary 
conditions ores 

+ (?."V) 8\/9xay^ » 0 

■» o 

3.1.2 Solution Methods 

The WECAN program can solvt equations 1 to 7 for any arbitrary 
temperature distributlcrM (The problem of designing a furnace configu- 
ration to produce a desired temperature distribution will be considered 
after it is learned what the desired distribution is.) Unfortunately, 
finding a solution for equations 1 to 3 does not solve ths buckling 
problem. In fact, the trivial solution, w s 0, is such a solution; it 
represents the nonbuckled suate. As long as no other solution exists, 
there can be no buckling. The problem thus is to determine the range of 
temperature distributions which induce sufficiently low stresses such 
that equations 3 to 7 have a unique solution, w =* 0. This range of 
temperature profiles will give leeway for proper silicon growth. We 
shall discuss four possible ways to approach this problerai; 

1. Tl\e WECAN computer program has the capability to solve 

eigenvalue problems, and we can approximate the buckling problem 
as an eigenvalue problem by varying the stresses in a fixed 
proportion; l.e., let 


a » 

\a o 

XX 

XX 

a » 

Xo° 

yy 

yy 

0 “ 

Xoo 

xy 

xy 


where ^xx, 'vy* ^xy are the stresses determined from equations 1 and 
2 for some typical temperature distribution T° (x) . WECAN can then 


8 


ORIGINAL PAGI \k 
OH POOR QUAUTV 


detemino thcs vuluos of A fotr which equations 3 to 7 have n unique 
trivial aolutlen. For those values of A, any of the toniporaturo 
dlsurlbutlona AT®(x) + ax 'h h (a and b are Independent of x and 
chosen to satisfy appropclato boundary conditions), would provide 
bucklo“freo silicon web growth. Notice that the oecond derivative 
of any of these distributions is proportional to that of the initial 
temperature distribution. Repetition of the above method for one or 
more initial temperature distributions (which are not proportional) 
could give some indication of the range of temperature distributions 
for which web growth would be stable. 

2. Timoshenko and Gere^^^ show how the eigenvalue problem mentioned 
in equation 1 can be solved as a minimization problem when 
aeries solutions to equations 3 to 7 can bo determined. They 
give sorleo expansions for similar boundary problems, and it is 
possible that these might be adapted to the problem at hand. 

3. Similar to the energy method Implied in equation 2 Is functional 
analysis theory^^'^ for the uniqueness of solutions to the 
buckling different tnl equation. This theory may provide 
additional Insight for determining the range of desirable 
temperature distributions for web growth. 

A. Integral equations provide powerful methods for obtaining 

solutions to equation 3^^^; however, the complicated boundary 
conditions, equations A to 7, make the analysis too difficult to 
pursue for the present. 

Conclusion ; Although the analytic techniques investigated offer 
great power and generality for establishing buckling criteria, they 
would require considerable additional development before they could be 
applied to our specific problem. Thus, for the present, it would appear 
that finite element numerical analysis, while less general in scope, is 
nevertheless the most effective approach for investigating thermal 
stress and critical buckling conditions in dendritic web crystals. 
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3.1.3 8 troB 0 ami lluckling Calnilatlono 


3. 1.3.1 WRGAN Ci)mpuloi‘ Codo 

The WeotlngUouBO finite element eodo, WlCAll, Imo boon need 
pcovlouBly for the rolntlvoly nlmplo problem of cnlenlntlng thermal 
sferoaaeQ in o two-dimenBlonnl reproaontntion of dendritic wob» In the 
proBont work> wo are employing the full copnbllitlea of the code to 
calculate the critical buckling atreoBes in a web crystal modeled as a 
throe“dlmonoionol bar including the stiffening effect of the supporting 
dendrites. The required three“dlmonslonal program uses nodal points at 
the intersection points of each of the finite element parallelepipeds 
and also at points midway between them. F\irthor, for the buckling 
analysis each nodal point is allowed three degrees of freedom. 

TO reduce computer time and cose, the grid else can be minimised 
by dividing the ribbon into only the number of finite elements needed 
for roasonable computational accuracy. While this iwmber cannot bo 
determined a priori, at least the plneement of the finite elements 
should ho arranged so ns to have the finest grid structure in the region 
of maximum temperature and s tress variation* For the dendritic web 
problem, this region is near the melt surface. A nonuniform grid can 
make the variation of stress across each finite element more nearly «ero 
tl»an can bo done with a uniform grid. Although WFXAN cannot handle 
large variations in grid spacing from one element to the next, it can 
accommodate a grid with element slse varying ns a geometric progression. 
Such a grid Cor the silicon web is shown in Figure 1; the elements next 
to the aurfaco are one-eighth as long as the elements at the other 
end. A perspective view of a section of the mesh la Illustrated iu 
Figure 2 to show tlve modeling of the bounding doudrltes. 

In addition to the finite element grid, WKGAN requires a 
temperature input for each nodal point. The temperature distribution 
along the ribbon cannot be measured easily but is obtained by into” 
grating the heat coviductlon equation for a specific furnace lid and 
shield geometry. The uniform step sice vised in the previous 
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Pigure 1 — Grid Used to Model Silicon Web for Stress Analysis 





Ifttcgratlon algorithm dooa not match with the geometrically varying grid 
uned with WBCAN, and modification to match the points led to numerical 
instabilities. A solution was finally obtained by Interpolation of the 
results of the uniform step integration. 

3. 1.3. 2 Buckling Analysis 

Buckling occurs as a sudden deformation in the planar silicon 
strip when the thermal stresses reach some critical mognitude. Mathe- 
matically, buckling can be represented as an eigenvalue problem where 
the critical stress distribution is represented by an eigenvalue and the 
buckled shape by an eigenvector, for buckling problems, WEGAN t^nlves an 
equation of the form 


[K] (A) « \ [Kg] {A} 


(d) 


wtiere [K] is the conventional stiffness matrix, [Kg] is the initial 
stress matrix resulting from a given temperature distribution, and (A} 
is a vector containing the unknown displacements of the ribbon surface. 

A solution parameter to the above equation, A, la nn eigenvalue and the 
corresponding (A) is an eigenvector; there can be many eigenvalues for n 
given equation. The smollesU positive one is the most significant, 
since the product of it with [Kg] is the weakest stress distribution for 
which buckling occurs. Negative eigenvalues imply a reversal of the 
sign of the stress distribution or a heating Instead of a cooling of the 
ribbon from the melt surface. Although negative eigenvalues may be 
obtained mathematically for some [Kq 1» they are not physically 
meaningful. We shall discuss these concepts further in terms of a 
specific example — ribbon grown using the J181 lid and shield 
configuration. 
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3 *1*3. 3 Sfci '080 and Buckling In J18X Ribbon 

J181 cibbon la gvown oufc of a nllleon welt covered and aUielded 
an illustrated in Plguro 3. The corresponding themial wodel Is shown in 
Figure the temperatures are measured datoi^^^ A previously developed 
thermal modeling program^^*^^ was used to calculate the toraporature 
distribution In the growing web resulting from this growth geometry} the 
result Is shown In Figure 5. This temperature distribution was used os 

the Input to WRGAN to calculate the stresses In a 10 cm long and 150 nn 

thick ribbon. The width of the web between the dendrites was taken as 
2.5 cm for a total ribbon width of 2.7 cm. The WECAN code then calcu- 
latod the x stress along the length of the ribbon, the y stress across 

the width of the ribbon, and the xy shear stress. Figure 6 shows the 

stress contours Cor these three components. The top edge of the strip 
Is tke outer edge of the dendrite, the left edge is the molt surface, 
the right edge is 10 cri from the melt surface, and the bottom edge Is 
the centerline of the ribbon with the s tresses symmetric about this 
line. Those stresses determine the matrix [Kg] of equation 1, since 
other components of the stress tenor vanish Cor thin ribbons. 

Using the initial sfeeos matrix [KgJ, WISGAN then solved equation 
13 Cor the eigenvalue X » 1.65, and li s corresponding eigenvector or 
buckled shape is illustrated in Figure 7. Since X is greater than 
unity, this ribbon would not have buckled all the stresses would have 
to be increased by a factor of 1.65 to obtain the illustrated deforma- 
tion. The stability of 150 pii thick and 27 mm wide web agrees with 
growth experience with the J181 configuration. 

In the actual growth apparatus, the ribbon is actually about 50 
cm long, rather than 10 cm, before It is wound up on the reel. The 


Only temperature for distances less than 10 cm from the melt was 
calculated using the previous program. An asymptotic procedure, to 
be discussed later, was used for greater distances. 
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FlTjure 6 — Stress Distribution in 27 na "J181 
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Figure 7 — Calculated Shape of Buckled Web (half strip) * 
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eomputer uoed foe the WEGAN analyolo dooo not have oufficient otoeago 
eapncity to model thio length> unloos It weeo divided up into very large 
eloraonts. Sln'^o the laegoot oteoeoeo occur within only about 3 cm of 
the melt interface, it was felt that a 10 era long model would be 
adequate for the present purposes. Nevertheless, wo added throe rows of 
cleroonts for a total length of 13,6 cm to check the invariance of the 
eigenvalue 1. Tliis longer ribbon was found to require slightly more 
stress for buckling with X ■ 1,67. We felt that thio 0.02 increase in X 
was within experimental error and insufficient to warrant the additional 
computer time. 


Not only is more computer time necessary to evaluate longer 
models, but the temperature input dota is more difficult to obtain. The 
temperature integration becomes inaccurate at large distances away from 


Che melt surface. An asjnaptotic racthed was used to find 
of the ribbon away from the molt, as graphed in Figure 5. 
wlil be discussed in an appendix to a future report* 


LUU vciiiijCLa 

This method 


3. 1.3. 4 Investigation of Buckling Porameters 

With some assurance that the WECAN buckling model is reasonable, 
we proceeded to use it to investigate the effect of ribbon width and 
thickness on the buckling resistance. These parameters are of prime 
Iraportance in the eventual goal of designing a growth geometry which 
will grow wide, thin (l»e., fast) web crystals. 

Width Effects . We chose two widths, 33 mm and 39.5 mm, to 
compare with the 27 mm J181 ribbon which was initially modeled. The 
stress contours for these widths are shown in Figures 8 and 9. Accord- 
ing to Boley and Wiener’s approximation, the shear stress should 
vanish on the centerline and it does in the numerical results. The x~ 
and y-stresses, however, do not exhibit the expected width dependence? 

and Oy~ w^ . Examining the point of maximum stress along the 
centerline, we find Instead that and Oy'^ w^*^® with a slight 

variation in the exponent depending on which of the width cases are 
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Figure 9 — Stress Distributions for 39.5 rm "J181" Web 




compared. The fact that the exponent fa not equal to 2 for the x-atreaa 
la not too Burprlsing oinco Holey and Wletior'a approxlmation^^^^ is best 
at diatanceo far away from the ends of a ribbon, whsreaa in our caae the 
maximum etroaa la only about a ribbon width away from the molt end. For 
this same reaaon, it cannot be expected that the x-streas at the den- 
drite edso ahould be double that at the center] ine and oppoaite in sign 
as predicted by the approximate theory. It is also true that the 
dendrites affect the stress pattern. The ratios of the minimum x-stress 
on the dendrite edge to the maximum x-stress on the centerline are -1.7, 
-1.8, and -1.9 for the 27 mm, 33 nun, and 39.5 mm widths of ribbon, 
respectively. The increasing magnitude of this ratio for increasing 
widths probably reflects the decreasing stiffness effect of the dendrite 
since the dendrite width is held constant at I ram* 

The smallest positive eigenvalues for the 33 and 39.5 mm 
ribbons are X « 0.797 and X » 0.A18, respectively. Their corresponding 
eigenvectors (shapes) are very similar to that shorn in Figure 7. Since 
these eigenvalues are less than unity, the ribbon would buckle before 
Its width could be increased from 27 mm to 33 mm. If we assume that X 
varies inversely as the Uj.j^ power of the width, then n “ 3.63 for a 
comparison of the 27 mm and 33 mm widths; using n « 3.63 to Interpolate, 
we find tliat the ribbon should buckle (X “ 1) at w =■ 31 min. This agrees 
well with growth experience for web crystals grown from a Jl81 
configuration. 

The 39.5 mm ribbon has another eigenvalue less than unity: X “ 

0.95. Its buckled shape is pictured in Figure 10. This wide ribbon 
could buckle like this or in the shape of its first order mode (X » 
0.418). The first order mode is the more likely one to occur with the 
higher order mode occurlng only if physical constraints should damp the 
lower mode. 

Thickness Effects . We evaluated two additional ribbon slaes; 

100 pm thick by 27 mm wide and 300 pn thick by 39.5 mm wide. From plots 
of the stress distributions, we find that the ratios of minimum to 
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Figure 10 — Buckled Shape for 39.5 ma "J181” Web (second eigenvalue) 


maximum x-atcesa are "1.4 for the araallor and -2.3 for the larger 
ribbon. In the 300 jm by 39.5 mm ribbon, the minimum x-atreaa does not 
occur at the anmo distance from the melt surface as that of the maximum 
x-8treso« For this cose, if wo choose a point on the dendrite at the 
same distance from the melt as the point where the maximum x-stress 
occurs on the centerline, then the ratio of dendrite to centerline 
stress is -2.1, which is very close to the predicted value of -2. Again 
we find that the ratio of dendrite to centerline stress is closest to -2 
in those ribbons in which the dendrites have the least volume relative 
to the rest of the ribbon. The difference between tlie thickness of the 
dendrite and that of the ribbon remains constant as tlie ribbon thickness 
is varied, so that the dendrite’s relative volume is smallest for wide, 
thick ribbons. 

The smallest eigenvalue for the 100 pm by 27 imn ribbon is l.ll, 
while that for the 300 pn by 39.5 mm ribbon is 1.06. Thus, the 27 mm 
wide ribbon would buckle when its thickness decreases to some value less 
than 100 pn, and the 39.5 mm ribbon would buckle at some critical value 
between 150 pm and 300 pn. To determine the critical thickness, wo 
interpolate using the assumption that the eigenvalue \ is proportional 
to the mth power of the thickness; for our cases, m =■ 0.98 for the 27 mm 
ribbon and m “ 1.34 for the 39.5 mm wide ribbon. Using these m’s, we 
find that the 27 mm ribbon should buckle at 90 pm and the 39.5 mm ribbon 
should buckle at 287 pn. These points, together v^ith the 150 pn 
thickness for a 31 mm width are plotted in Figure 11. The regions of 
stability (thick or narrow web) and instability (wide or thin web) are 
separated by an essentially straight line. Points determined by 
calculating the width at which a 100 pt or 300 pn thick web would buckle 
also lie on this line. The simplicity of this result is rather 
surprising in light of the complexity of the buckling mechanism. 

3,1.4 Generalized Design Guides 

Tlie preceding sections have presented the development and 
application of finite element methods for calculating thermal stresses 
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Figure 11 ““ Buckling as Function of Width and Thickness for J181 Lid 
Configuration 
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and also critical buckling conditions during the growth of silicon 
dendritic web. Evaluation of the buckling behavior of a given growth 
configurotion is a three-step process! (1) calculation of the 
temperature profile in the web crystal using the temperature thermal 
code; (2) calculation of thermal stresses for a three-dimensional model; 
and (3) calculation of the buckling eigenvalues. In practice, steps 2 
and 3 are performed in the some computer program; however, this computa- 
tion is rather lengthy and complex. By comparison, calculating the 
stresses in a two-dimensional model is relatively simple. It would be 
desirable to have some qualitative or semlquantitative guidelines to 
evaluate a growth configuration based only on two-dimensional stress 
calculations. This section presents some preliminary work on such 
guidelines both for residual stress and for critical buckling stresses. 
The results are illustrated by consideration of the J98M3/J352 
configuration. 


3. 1.4.1 Residual Stress Index 

Up to now, our analysis has emphasized thermal stress buckling 
due to its limiting effect on ribbon width. Residual stress, on the 
other hand, may not as directly affect the ribbon growth but it con 
certainly affect the ease of fabrication of the ribbon into devices. 
Consequently, we must also consider residual stress when evaluating 
various lid and shield configurations. 

Although the generation of residual stress during crystal growth 
is a complicated viaco-elastlc process, hopefully a simple model will be 
sufficient to predict growth behavior, at least on a comparative basis. 
In simple terms, residual stress occurs when the thermal stress exceeds 
the yield stress for silicon. In the hlgh-teroperature region near the 
growth front, the transverse stress, Oy, is the only stress component of 
any magnitude, and we shall use the following empirical relation for the 
yield stress^^^^: 
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o^g 2.57 X exp (49459/T) dytt/ci«2 


( 9 ) 


where T is the teraporature in Kelvin. 

The amount of residual stroBa should depend on the maenltude of 
the difference between the lateral streaa and the yield otresa and also 
on the length of time that the material ia in that condition. Thus, to 
obtain some approximation to the residual stress, we define a Residual 
Stress Index (RSI) as the area between the yield and transverse stroaa 
curves plotted as functions of distance from the melt. For this 
purpose, the transverse atreas at the centerline is used since it is 
largest in this region. Figure 12 illustrates the area concerned using 
case 10**6 (next section) as an example. 

For computational purposes, the stress functions are approxi- 
mated by line segments connecting data points; (xi,yi) and (xi,a£)i 1 ^ 
0, ..., n. The are distances from the melt, the yi's are the 

ordinates of Che yield stress, and the z^’s ere Che ordinates of the 
transverse stress; n is the smallest value of i such that lyjLl>Zil. The 
Residual Stress Index may thus be computed as the sum of the areas of 
trapezoids and a triangle. 


n-1 


RSI = 1/2 (Kj - x^_^) 


+ 1/2 (y„_i / (yn-r^n + 


In the case illustrated, the RSI “ 41.5 Mdyn/cmj the magnitude 
of this nun\ber indicates the likelihood of undesirable residual stress. 

3. 1.4. 2 Stress and Buckling 

It would be convenient to have an index for buckling analogous 
to Che Residual Stress Index. No single number, however, could be 

/■*% 
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Figure 12 — Schematic of Integration for Residual Stress Index (RSI) 
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expected to predict buckling in general* 1110 otrooo patterns determine 
the buckling for a given ribbon alise. With the goal of understanding 
which foatureo of the atreao pattorna are moat Important, we examined 
the variational formula for the buckling eigenvalue, 

P / J " yy>^ - 2(1 - U) - »^y^) 1 tlxd, 

/ / + Ny«/ + dxdy 


where D la the flexural rigidity of the ribbon, \> la Poisson's ratio, 

N^, Ny, and N^y are the longitudinal, transverse, and shear atreaaea, 
respectively, which have been multiplied by the thlckneaa of the ribbon, 
and S “ w(x,y) la the buckled shape of the ribbon or, equivalently, the 
surface which aatiafleo the appropriate boundary conditions and minimizes 
the expression 11* Tl\c subscripts on w refer to its pvirtial derivatives 
with respect to the indicated variables. 

We examined the values of the denominator Cor some of the 
buckling models of the J181 growth configuration and found that the terra 
Involving the shear stress, ^N^yW^^Wy, was the major contribution to the 
Integral and that this term was greatest near the melt surface. 

Although the terra Involving the x and y stresses directly are small (but 
not negligible), these stresses Influence the buckled shape w(xy), and 
thus affect the term Involving the shear stress. 

As an example, compare the stress data in Figure 6 (Case lA) 
with the data in Figure 13 (Case 2-2*50). Case lA presents the results 
for a .015 cm thick ribbon, while Case 2-2.50 represents a .010 era thick 
rlbbonj otherwise the two cases are identical. If the stress patterns 
were Identical, then equation 11 would requite the buckling eigenvalues 
to be proportional to the square of the ribbon thickness} D has a factor 
of the cube of the ribbon thickness, while the N's have a single factor 
of the thickness. Thus, the eigenvalue of the .015 cm ribbon (Case lAt 
X * 1.65) should be 2.25 tiroes as large as the eigenvalue for the .010 cm 
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ribbon (Case 2-2.50; X - l.ll). Tlve finite clement buckling calcula- 
tions ohow that this ratio is only 1.49. 


One reason that the theoretical ratio of 2.25 is probably too 
large is that the dendrite size has been kept the some in both cases 
(.06 cm thick) and tends to dominate the stiffness factor, D. Even 
taking an average thickness gives a ratio of 1.79, still too large. 
Furthermore, the stress patterns are not Identical for both cases, the 
temperature profile is slightly different for the two ribbon thicknesses 
and causes a slight change in the stress distributions. The differences 
in the maximum to minimum x-stress, y-stress, and xy-stress has changed 
by 8X, -4%, and -1%, respectively, for the .015 cm case, compared with 
the .010 cm case. The apparently most significant factor In the denomi- 


-t ft 


tivs shear stress, Increased in 


nacof dr equatiuu iii, the large nega 
magnitude by only Z%» This alone does not account for the ratio of the 
two eigenvalues; the changes in the stress patterns have changed the 
buckled shape significantly. Since the x-stress pattern had shown the 
largest increase In peaks and valleys, it is the most likely contributor 
to the larger-than-expected potential for buckling (smaller eigenvalue). 
The maximum and minimum x-s tresses occur at the ribbon center line and 
at the dendrite and at about the same distance from the melt. Thus, the 
difference between the maximum and minimum x-stress is almost the same 
as the maximum of the Ax-stress defined as the difference between the x- 
8 tress at the centerline of the ribbon and at the dendrite at the same 


distance from the melt. 


For the J352/J98M3C lld/shield configuration which will now be 
considered, the Ax-stress, x-stress, xy-stress, and the y-stress (away 
from the interface) patterns show approximately the same behavior with 
respect to the parameter variations under consideration* In other 
words, even though the Ax-stress pattern is not the only contribution to 
buckling, It is a good indicator for it since all the other stress 
patterns increase or decrease along with It for the cases we examined. 
While the y-stress near the melt is important for residual stress, it 
appears to be much less so for buckling. 
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3tl.4t3 Application to the J352/J98M3C Growth ConfigurationB 

The J352 and J98M3G growth configurations are two lld/shleld 
configurations which have been used a number of times to grow silicon 
web crystals. Both configurations have some useful properties for web 
growth! the J352 produces crystals having very low residual stress, 
while the J98M3C, a width-limiting lid, produces constant width crystals 
having moderately small residual stress. 

The web temperature profile is calculated as a one-dlmonaional 
problem, so that only the vertical geometry of the lid and shields Is 
important. Furthermore, the resolution of our present thermal model, 
shown in Figure 14, Is such that both the J352 and J98M3C configurations 
can be modeled by more or less the same geometry. In the model, the lid 
itself is more or less accurately represented as to thickness, slot 
width, and possible bevel and may also have a different slot tempera- 
ture, TC, and top temperature, TS, a condition which would exist with a 
thin top shield very close to the top of the lid. In the model, 
however, the stack of shields above the lid must be represented by a 
lumped region of constant temperature, TAS, which must be chosen to 
approximate the effective temperature of the shield assembly. In the 
calculations which are discussed here, the geometrical factors were 
based on the actual lid and shield dimensions, while the effective 
temperatures were estimated from measurements of various lids and 
shields . 

Table 1 shows the values of the various model parameters for six 
cases related to the J352/J98M3C growth configurations. In addition to 
the physical parameters illustrated in Figure 14, there is a model 
parameter, V, which estimates the growth velocity of the web (the moving 
frame velocity in the heat conduction equation). This velocity depends 
not only on the heat lost from the web itself, but also the latent heat 
lost to the supercooled melt, and must therefore be estimated before the 
model can be run. One of the model output parameters, VV, is the 
portion of V due to the loss of latent heat through the web itself. In 
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Figure 14 ““ Schematic J35Z/J98M3C Configuration and Equivalent Thermal 
Model Parameters 
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addition to VV, the table also gives the Residual Stress Index, RSI, and 
the maxiraura Ax-a tress* 

Each of the cases listed in Table 1 hos certain specific 
features which cm bo discussed individually. 

Rerooval of top shield was simulated accidently in Case 10~l by a 
numerical error in ascribing on input value of W/^ *•7.9 cm instead of 
0.79 cm. This case modeled moving the top shield far from the lid, 
effectively removing it from the calculation. Case 10-3 is identical to 
Case 10-1 with the proper » .79. Comparison of the results shows 
several effects; 1) Stresses are larger and occur closer to the melt 
when the shield is removed, so that under that condition there would be 
an Increased tendency to buckle} 2) the residual stress index la larger, 
so that there would be greater plastic flow near the interface} and 
3) the web growth velocity is faster. 

Melt level changes are modeled in Cases 10-2 and 1-3} in Case 
10-2 the melt level is such that the growth front is at the same level 
as the bottom of the lid, while in Case 10-3 the melt level has dropped 
0.1 cm. A comparison of the model output parameters shows that for the 
lower melt, the web component of the growth velocity decreased by about 
20% and the maximum Ax-stress increased slightly, but that the Residual 
Stress Index changed by a factor of about 3.5. Thus, changes in melt 
height should have negligible effect on the buckling of the ribbon, a 
moderate effect on growth velocity (which has been observed) and a very 
pronounced effect on the residual stress. 

Sensitivity to assumed growth velocity was tested by Cases 10-2 
and 10-4. As discussed earlier, the moving frame velocity V Is required 
as an input parameter for the integration, whereas the web velocity W 
(“ K/pL dt/dx]jj„Q where k » thermal conductivity, p « density, and L “ 
latent heat of fusion) is calculated from the program. In Case 10-2, V 
was set at 2.0 cm/rain while in Case 10-4, V was set at 1.5 cra/min} 
comparison of the output parameters shows that there are only negligible 
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Table 1 — Model Parameters for J352/J98M3C 


Parameter 

lO-l 

10-2 

CASES 

10-3 

10-4 

10-5, 

10-6 

(Input) 

TS 

1450 

1450 

1450 

1450 

1450 

1575 

TC 

“K 

1625 

1625 

1625 

1625 

1625 

1625 

TA 

OR 

300 

300 

300 

300 

300 

300 

TAS 

“K 

1175 

1175 

1175 

1175 

1175 

1175 

TA2 


1450 

1450 

1450 

1450 

1450 

1450 

B 

era 

.015 

.015 

.015 

.015 

.015 

.015 

W1 

era 

.3 

.3 

.3 

.3 

.3 

.3 

W2 

era 

3.1 

3.1 

3.1 

3.1 

3.1 

3.1 

W3 

era 

4.5 

4.5 

4.5 

4.5 

4.5 

4.5 

W4 

era 

7.9 

.79 

.79 

.79 

.79 

.79 

W6 

era 

4.5 

4.5 

4.5 

4.5 

4.5 

4.5 

ZD 

era 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

ZH 

era 

1.75 

1.75 

1.75 

1.75 

1.75 

1.75 

TLID 

cm 

.8 

.8 

.8 

.8 

.8 

.46 

TBEVEL 

cm 

.25 

.25 

.25 

.25 

0 

0 

LIN 

cm 

.1 

0 

.1 

0 

0 

0 

V 

cm/min 

2.0 

2.0 

2.0 

1.5 

2.0 

2.0 

(Output) 

W cra/min 

1.20 

1.39 

1.12 

1.40 

1,09 

1.68 

RSI 


33.9 

12.4 

3.5 

13.1 

3,2 

41.5 

Max ha 

X 

Mdyn/cra^ 

782 

645 

654 

665 

716 

528 
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changes; hence, for velocltlos currently being calculated, the 
complication of Iterative solutions Is not warranted. 

Lid bevel effects were brought out by comparison of case 10-5 
wl^’h the standard case 10-2. Both geometries were Identical except that 
the lid In Case 10-5 did not have the .25 cm x 45® bevel of Cose 10-2. 

The effect of having the otralglit-alded slot was to slow down the web 

velocity, VV, and also to decrease markedly the Residual Stress Index. 
There was also some increase in the maximum Ax-stress, but this appears 
to be relatively small compared to the other changes. 

Lid thickness was the final parameter Investigated. Case 10-6 

had the same geometry as Case 10-5, except that the lid was thinner; 

0*46 cm instead of 0.8 cm. The result was that the thinner lid showed a 
marked increase In VV and also a large increase in the Residual Stress 
Index. The maximum x-s tress decreased but, in addition to the main 
peak, developed a secondary peak nearer to the growth front. This 
second peak complicates interpretation of the results In terms of the 
tendency to buckle, and a complete buckling calculation would be 
necessary to evaluate the change. 

Six conclusions can be stated from tlie modeling runs for the 
J352/J98M3C model reported in the previous section: 

1 . The presence of a top shield apparently reduces the thermal 
stresses in the ribbon — both the RSI and tlie Ax-stress. 

2. Lower malt level gives a lower RSI but does not affect the x- 
stresses further up the ribbon. 

3. The calculations are insensitive to the pull velocity assumed as 
an input parameter, at least for a reasonable range of values. 

4. The presence of a bevel on the susceptor lid can have a signifi- 
cant effect on the partial web velocity, VV, and on the thermal 
stresses . 

5. A thin lid Is likely to have a large RSI, albeit a larger web 
velocity. 
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6* The growth velocity and the thermal etreoaea are unrelated, 

although they may appear to be correlated if only one parameter 
is changed. 

3.1.5 Zero Stress Profiles 

We have seen that the J181 configuration produces a temperature 
profile along the ribbon that would cause sufficient stress to buckle 
silicon ribbon grown to a 31 mm width and 150 pm thickness. In order to 
redesign the lid and shield configuration to produce wider and/or 
thinner ribbon, we would like to know how the temperature profile could 
be modified to reduce stresses in the ribbon. A first step toward this 
goal is to find temperature profiles that produce no stress at all. 

That a linear temperature profile should yield zero stress is 
derived from the equation 


a 2<|> « -a EAT 


( 2 ) 


where is Airy's stress function, a is the thermal expansion coeffi- 
eient, E is Young's modulus, and A is Lapalace's operator. For zero 
stress, the left-hand side vanishes and we find that T"(x) ” 0, since T 
is a function only of the distance x from the melt and not of y or z (to 
a close approximation in ribbon growth). The above equation, however, 
was derived for constant a. For a a function of T (« => 2.8171 x 10“^ 

+ 9.789 X IO^IOt for silicon), we replace the condition T" « 0 by (a T)" 
» 0. One resulting temperature profile is shown in Figure 15, curve 9- 
1. This curve must pass through the point (x » o, T “ IbSS^K) and can 
be at any chosen temperature at the other end. However, if this 
temperature is reasonable, then the slope at x « 0 would be. so small as 
to produce very slow ribbon growth. If the initial slope of this curve 
was chosen for reasonable speed, then the tempefiiture at the otlier end 
of the ribbon would be so low as to be physically impossible. One way 
around this problem is to try to design for a "zero stress" temperature 
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profile near the molt with a reaoonablo slope and then, at ooroe point 
far from the melt, decrease the temperature at a much slower rate than 
that of the "zero stress" profile. 

In Figure 15, curve 9-1 is a "zero stress" curve. At 8 cm from 
the melt, two "toils" were joined to this curve in such a way as to 
maintain a continuous first derivative. The "tail" 9»*3 is linear and 9- 
2 is an exponential function of x. 

Plots of the difference in longitudinal stress between the 

centerline and the dendrite edge for each of the three temperature 
profiles are shown in Figure 16. The linear tail, which is barely 
distinguishable from the "zero stress" tail, produces a small additional 
stress. The exponential tall produces a major stress peak at about 8.5 
cm from the melt, although it is not nearly as strong as some of the 
J181 cases which peaked at close to 1000 Mdyn/cm^. The "zero stress" 
profile appears to produce small stresses at both ends of the rlbbonj 
this is an artifact of the model. For these three cases, an 8 x 100 
mesh of linear, two-dimensional, finite elements was used to model the 
half-strip. The "zero stress" case was additionally modeled with a 
5 X 23 mesh of cubic finite elements which have eight more nodes per 
element than the linear elements. The calculations using the cubic 
elem-vints gave a zero stress over the whole ribbon length and Incldently 
took much less computer time. 

We have found that small deviations from the ideal "zero stress" 
temperature profile can produce major stresses in a ribbon. To under- 
stand and quantify this effect better, we plan to examine additional 
synthetic temperature profiles for stress generation. In particular, 
profiles with steeper slopes at the melt surface are required for 
reasonable growth speeds. 

3.1.6 Temperature Data for Modeling Calculations 

Thermal modeling of a web growth system requires lid and shield 
temperatures as input parameters to calculate the heat transfer between 
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Figure 16 — Comparison of Stress for "Ideal" and "Near Ideal" 
Temperature Profiles 
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the web and Its surroundings# In most modeling work done so far, the 
necessary temperature data has been estimated from a limited set of 
measurements made on a J18X configuration some time ago.^^^ Current 
modeling has become increasingly sophisticated so that more accurate and 
complete temperature data are needed. Not only are more extensive data 
needed, but also information is required as to the variation in 
temperature which can be obtained by adjusting the coil height and 
changing the coupling of the shields to the work coll, etc. 

Two current growth configurations, the J352 and the J98M3A, were 
chosen for measurement since they are representative of many other 
configurations. The general lid and shield geometries are shown in 
Figure 17, which Illustrates the configurations and shows the thermo- 
couple locations. Additionally, two types of shields were used as the 
top shields in the J352 configuration: a normal shield which had a 

continuous edge and a slitted shield that had 1 x 12 ram slits every 
12 ram along the periphery to reduce coupling with the work coll. 

Previous measurements had shown that such a slitted top shield ran 
significantly cooler than a similar solid shield, and additional 
information was needed. 

3. 1.6.1 Experimental Procedure 

Shield temperatures, which were generally less than 1300°C, were 
measured with Type-K (chromel-alumel) thermocouples, while the lid 
temperature was measured with a Type-B (Pt/30%Rh — Pt/6%Rh) couple. 
Small holes were drilled In both the lid and the shields and the tliermo- 
couple beads cemented in place with alumina cement to ensure good 
thermal contact. The thermocouple leads were brought through the 
furnace chamber wall with Conax compression seals. The Type-K leads 
were continuous to the measurement apparatus, while hlgh-purlty copper 
was used for extension leads with the nobel metal couple. The thermal 
emf’s were measured and converted to temperature data by a Flukve 2240B 
Data Logger which also printed out the results at regular intervals. 
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Figure 17 — Configurations for System Temperature Measurement 
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Moot of the teraperatuiCG data were taken without silicon in the 
system so tlmt n wide range of susceptor temperatures could be covered# 
For each configuration, temperature data wore taken at three different 
coil heights and a range of susceptor temperatures (measured by a I'ype-B 
thermocouple on the axis of the susceptor) . For the case of the J983A 
configuration, additional data were taken with a normol silicon charge 
in the system so that the lid and shield temperatures could be measured 
at the "hold” temperature to standardize the temperatures. 

3. 1.6. 2 Results and Analysis 

Measurements made without silicon in the system showed that the 
lid and shield temperatures wei*e related linearly to the susceptor 
temperature over a wide range. Generally, the shield temperatures did 
not change as rapidly os the susceptor temperature; the lid temperature 
had nearly the same rate of change as the susceptor. Table 2 summarizes 
the measurements made on the various configurations shown in Figure 1. 
The data are expressed in terras of an ampirlcal equation: 

T » A + B(Tg“1380) I°G), where T is the element temperature, Tg is the 
axial susceotor temperature, and A and B are the fitting coefficients. 
The form of the equation was chosen so that A represents an approxima- 
tion to the element temperature near normal operating conditions. 

The coil position obviously affects all the element temperatures 
but has the least Influence on the lid temperature. Slitting the edge 
of the top shield also has a strong effect on the shield temperatures, 
which shows that coupling with the work coll is a very Important factor 
In determining the system temperatures. We attempted to analyze the lid 
and shield temperature distribution in terms of radiative transfer with 
an effective eraissivity as the adjustable parameter. An emissivity on 
the order of 1.5 was needed to give even an approximate fit, again 
indicating that the inductive coupling was indeed a very Important 
factor, even with the slit ted shields. The results thus confirmed our 
assumption that measured, rather than calculated, temperatures were 
essential for accurate modeling predictions with real systems. 
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Table 2 


Summary cf System Temperature Data (Without Silicon) 
T - A + B (Tg - 1380) 


Bottom Shield 
3 A B 


1.0106 

.8968 


Top Shi eld 

rn nil i> II in I II III III 

A B 

“C 


1031.4 

990.1 

973.0 


0.7031 

.7432 

.3456 


Configuration 


J352, 2 shields 
(top solid) 


J352, 2 shields 
(top slitted) 


1352, 3 shields 
(top slitted) 


J98M3A, 4 shields 
(top solid) 


Coil Lid 

Height A 


ram* 

“C 

93.0 

1324.9 

97.5 

1313.3 

102,0 

1296.1 

93.0 

1320.1 

97.5 

1306.4 

102.0 

1293.6 

93.0 

1328,2 

97.5 

1319.5 

102.0 

1306.6 

97.5 

1334.7 

97.5 

1322.3 

102.0 

1306.5 


.7988 

.8790 

.8687 

.8918 


.9723 

1144.8 

.6573 

.8702 

1112.9 

,6456 

.8506 

1075.7 

.4183 

.8694 

1177.3 

.7348 

.9379 

1144.3 

.7551 

.9576 

1107.5 

.6202 


942.4 

.7568 

936.5 

.8727 

908.0 

.7448 

944.7 

.7146 

918.7 

.7209 

883.2 

.5618 

1033.7 

.6174 

990.1 

.5910 

936.3 

.6573 


higher coil 


*Coil height measured from arbitrary zero; smaller numbers are 
positilons with respect to susceptor. 

Tg is the axial susceptor temperature. 


IvxpftclraentB wlUiout silicon in the eystera give the data required 
to determine how the element teinporntureo change when the susceptor 
temperature changes. The presence of; silicon permits evaluation of the 
element temperatures for actual operating conditions. Table 3 lists the 
lid and shield temperature data for a normal silicon melt at "hold” 
temperature for both the J98M3A configuration evaluated in the present 
experiments as well as the previously reported data for the J181 
configuration (a three-shield configuration) The data clearly show 
that the lid temperature has a relatively small range while the shield 
temperatures can vary widely. For example, in the J98M3A data, a 9 mm 
variation in coll height changes the lid temperature by only 17.6^C, 
while the shield temperature changes by 92.3°C. Furthermore, thicker 
lids run somewhat hotter: the 9.5 mm lid of the J181 runs about 25°C 
hotter than the 6. A kks lid of the J98M3A. Thus, lid thlcknoss provides 
a temperature parameter as well os o geometrical parameter for future 
system design. 

3.2 Experimental Web Growth 
3.2.1 Advanced Growth Studies 

The objective of this activity Is to develop the experimental 
understanding of web growth behavior in order to grow web crystals wider 
and faster, with the aim of achieving area growth rates up to 35 cra^/mln. 
This effort is tied closely to the thermal stress modeling effort. 

There are two major aspects to this work; (1) developing lid 
and shield configurations giving higher growth speed at greater width 
and (2) developing the capability for steady state growth for long times. 

The first activity is concerned with the continued design and 
testing of growth configurations which will permit both high-speed and 
low-stress operation, particularly low elastic stress which presently 
limits width through deformation of wide crystals. The second activity 
not only involves the design and validation of automated growth equip- 
ment, but also Includes the achievement of growth configurations which 
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Table 3 


Sumwary of Syatem Temperature Data (With Silicon) 


Configuration 

Coil 

Height 

ram* 

Lid 

Temperature 

Top Shield 
Temperature 
“C 

J98M3A, 4 shields 
(top solid) 

93.0 

1342.4 

994,9 

97,5 

1332,7 

946.0 


ior,o 

1325.1 

902,6 

J181, 3 shields 

(top solid) 

100,0 

1354.7 

963,1 


102.7 

1351.3 

941,7 


105,0 

1344,6 

927.4 

(top slitted) 

101,0 

Open TC 

876.6 


*Coll heiuht measured form arbitrary zero; smaller numbers are 
higher coil positions with respect to the susceptor. 
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produce a conscant, controllable ribbon width compatible with melt 
reploniohmont. 


3. 2 *1*1 Growth Velocity 


During the early part of this period, experimental efforto were 
directed toward gaining a bettor understanding of the configurational 
parameters which affect growth velocity. 


In order to Increase growth velocity, a series of experiments 
with thinner lids was undertaken. A thin lid has the effect of 
increasing radiative losses from the web near the growth interface so 
that the heat of fusion is more rapidly dissipated, thus permitting more 
rapid growth for a given crystal thick!:.'f;Ss. Previous work had shown 
viloC o*iiip*/ wiij.irikj.11^ v.nc iiu, witj-jLi: iiuvxuj^ tt xuxgu uxx«Cc on growcn 
velocity, can have negative effects in terms of crystal quality and 
oxide deposition. These effects occur because a thin lid tends to run 
colder in the region of the growth slot. Thus, concurrent changes must 
also be made in the top shield configuration in order to provide a 
vertical temperature distribution which not only produces substantially 
higher growth velocities, but also maintains good growth behavior in 
terras of growth initiation, crystal quality and freedom from oxide 
interference. To accomplish this, a series of experiments was undertaken 
in which the lid was systematically thinned in steps while changes were 
made in the top shield configuration to maintain good quality growth. 


Initial experiments involved an lid reduced in 

thickness by 20%, with the RE-1 shield arrangement. As expected, growth 
velocity increased significantly as determined by thickness-velocity 
measurements. In terms of the relationship 


V 


c 


+ 


d 

7t 


which relates growth velocity (V) to web thickness (t), both coefficients 
c and d Increase when compared to the original thicker lid designs. 
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However, growth initiation was difficult ond stability and crystal 
quality poor. The addition of two more top shields to this configura- 
tion resulted in stable growth with some trodeoff in maximum growth 
velocity. 

Some other variations in the number and spacing of top shielding 
ualn.'X the same lid have been tested. Only small effects wore found in 
terms of growth velocity, but these experiments hove led to the develop- 
ment of several Important guidelines relating to lid ond shield 
configuration which, if not considered, may produce unwanted side 
effects, for exomple, sporadic floating ice and growth-inhibiting oxide 
accumulation in the slot region of the shields. The configurational 
parameters to be controlled are the lid temperature and the spacing 
between the various components of the lid-shield stacks. By maintaining 
the bottom of the lid at sufficiently high temperature and keeping 
appropriate top shield spacing, oxide and ice are prevented. 

The lid temperature can be controlled to the required degree by 
adequate radiation shielding and by proper positioning of the work coll. 
In addition, the lid bottom and susceptor top must mate well to insure 
maximum contact and thus maximum heat transfer between the two. ^-Jhen 
conditions are properly controlled in terms of shielding, lid-susceptor 
fits, etc., the lid bottom is totally free of oxide deposits, and ice 
from spalling does not occur. 

One of the important aspects of top-shield spacing is the 
separation between the lowest shield and the lid. If this spacing is 
too small, oxide is deposited on the shields and grows out over the slot 
region, inhibiting crystal growth. The minimum spacing needed to 
prevent such oxide deposition appears to be about 4.5 mm. This 
observation had previously been made relative to the J-18X lid and 
shield configuration, in wlilch the lowest top shield is thin and 
relatively cold. Recent experiments indicate that this spacing 
requirement is much more general and that the temperature per se is not 
the only factor to consider; the convective gas flow patterns are 
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critical in relation to growth-inhibiting oxide deposition on the cop 
shields . 

Another observation relative to shield spacing relates to the 
spacing between the individual shields in multiple shield stacks. A 
spacing of 1.5 mm is generally satisfactory. Expanding this spacing to 
3 mm can cause the deposition of a hard blue oxide on the web crystal 
(as opposed to the normal, fluffy, easily removed type). Again, in this 
case the convective gas flow patterns are believed to be significant. 

Even with improved shield spaclngs and lids which had reasonably 
fast growth speeds, growth was not always easy and crystal quality 
sometimes poor. In some cases, the crystals would degenerate after some 
tens of centimeters of growth. Split width measurements showed that 
some of the crystals had considerable residual stress, even at the 
relatively narrow width of 20 mm or bo. These factors indicated that 
the lid temperatures in the vicinity of the slot were too cool, even 
with multiple top shields. 

Another approach to achieving the proper thermal geometry was 
then tried: a thick (9.6 ram) lid with a recessed center portion around 

the growth slot. Such a configuration would have a higher slot tempera- 
ture since the higher thermal conductance of the thick lid would 
transport heat more readily and, further, the thicker lid would couple 
to a greater extent with the work coll. 

Several recessed lids with thicknesses in the slot region 
ranging from 5.6 to 3.2 mm have been Investigated. The speed of these 
configurations is very similar to the speed of the flat lids of equal 
thickness, with the distinct advantage that the growth behavior is 
better. With each lid configuration a number of shield configurations 
has been tried. Two general types of configuration are shown in Figure 
18. In the first configuration flat shields are used, but in the second 
configuration, one of the shields is made of thinner material and bent 
over to reduce radiative heating from the sides of the recess. An 
advantageous configuration uses the bent shield together with a top 
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Typical Cross Sections of Recessed 
Enhancement 


Lids for Growth Speed 
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shield which lo alifeted niround Che edge to i*odvt?e coupling. Thla 
eonvhinntioii ptovidod tclntlvoly Caat gcowth wltV an acceptably low 
roaldual attouo, although other conCiguratlona were nearly aa good. 
Although the shield conflguratlona and Che slot ahnpe arc very important 
in determining Che opeod of the ayatem and its atreoa behavior, one of 
the most important paramotera ooema Co be the welt height. Later In 
this report (Figure 31), data arcs given tor the change in web thieknoaa 
as a function of crystal weight (which is equivalent to molt height). 
Using the CQVivorsion that 20 gm of crystal ia equivalent to I mm of melt 
height, it can bo seen that aignlf leant apeed improvement reaulUa when 
growth ia maintained with high molts. Some of the conflguratlona we 
tested have ahown prowlso of auch behavior. 

The growth experimenta with thin lids, recessed lids, and with a 
variety of shield geometries have roemphaslaed the large mimbor of 
para?neters which influence the growth behavior of dendritic web. All of 
these factors interact in a complex way to determine the growth speed 
and thermal stress In tl\o growing ribbon. Tlie experiments so far 
Indicate that achieving any optimum configuration by empirical methods 
alone is difficult. It should bo noted that the reduction in buckling 
stress through the utiliisatlon of the thermal modeling results wUl not 
only benefit crystal width, but also growth velocity through the thick’" 
ness dependence; l.o., low stress permits high crystal quality to he 
maintained at lower thickness, and thickness and velocity bear an 
inverse relationship. 

3 .2. 1.2 Width Control 

for the continuous growth of long web crystals, It is necessary 
to limit or control the width at some desired value. In addition to 
targeting width to a specific size solar cell (reducing waste), the 
crystal width can ho maintained below that at whicti stress'^lncluced 
buckling might occur. Thus» very long crystals can be grown with the 
majority of the material usable for cell fabrication. 
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Configurations for width control were developed on an in-house 
program and utilized in this program where applicable. Data from on 
early crystal for which the width was controlled to within a fraction of 
a millimeter by passive lid design techniques is shown in Figure 19. 
Crystals have been grown to lengths over five meters at controlled 
width, without the benefit of melt replenishment. 

3. 2. 1.3 Melt Replenishment 

During web growth, a constant melt level is one of the require- 
ments for maintaining constant thermal conditions, in turn required for 
steady state growth. Under the previous contract^ a method was 
developed for detecting the silicon molt level and maintaining it during 
web growth; the feasibility of controlled melt replenishment over 
periods of many hours was demons trated. On the basis of this experi- 
ence, it was clear that a number of refinements was desirable in order 
to improve the melt level control function and to meet the requirements 
better for long-term semi-automatic growth, a goal of the overall design 
of the new experimental web growth machine. This improved system is 
described in Section 3.3. In this section, we will briefly discuss some 
of the experimental work related to melt replenishment. 

During the period when the new melt control circuitry was being 
designed and constructed, growth experiments were carried out on the RE 
furnace using manually set replenishment rates. 

Lid and shield designs which provide width control were adapted 
for replenishment by the addition of feed and laser holes. This enabled 
the melt level to be maintained reasonably constant using a constant 
feed rate. Any mismatch between feed rate and consumption rate Is 
manifested as a thickness variation in the crystal, so that any needed 
adjustments in feed rate can be made empirically. Tlius, growth at 
constant vridth allows growth experiments to be carried out at nearly 
constant melt level without the need for melt level sensing, at least 
over periods of several hours. It was also established that the 
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Figure 19 — Variation in Web Width with Length of Crystal Growth for 
Lid Geometry Designed for Width Control 


width-limiting lid and ahlold doslgna contlnuod to carry out thoir 
function in a melt ropXcnlahment configuration. Crystala up to four 
meters long wore grown with manual replenishment. 

Although the standard rectangular crucible design has performed 
well, It has become Increasingly apparent that a longer crucible would 
be desirable In order to provide more physical space in which to 
separate the thermal requirements for growth of wide crystals and 
replenishment at the higher rates required to maintain n constant 
throughput. A new rounded-end crucible was designed and fabricated. 

This crucible has the same shape as the susceptor, l.c.) the geometry 
follows the susceptor perimeter. This design adds nearly 5 cm of useful 
length along the crucible center line with no change in t1\e furnace 
configuration except the shape of the recessed portion of the 
susceptor. Initial tests of growth behavior in the elongated crucible 
did not show any inherent problems, but the thermal charucterlatlcs of 
this configuration are somewhat different from those of the shorter 
rectangular configuration. Tlie required shield positions, coil height, 
and crucible barrier positions for good growth with slmultanoo\is 
replenishment were determined experimentally. The expected improvements 
in replenishment capabilities over the shorter crucible were indeed 
realized, l.e., Increased thermal Isolation of the food and growth 
regions of the crucible. Pellets melted easily at more than twice the 
feed rate possible in the short crucible configuration without affecting 
crystal growth. 

While this work was in progress, the new melt level control 
system was installed on the RE furnace and tested under growth condi- 
tions. Overall, ct\e system operated as designed. This was established 
by recording both the output from the detector amplifier and the output 
from the motor controller (input to the pellet feed motor) and observing 
changes with changes in melt level. One modification deemed desirable 
was to lengthen the time constant in the control system in order to 
reduce the noise generated by the shimmering melt. In order to provide 
sotae flexibility in the choice of time constant, a circuit was designed 
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which provldcB a eJioice of 0, 20, 100, or 400 boc. Tliio circuit was 
included in the welt level control system installed on the new experi- 
mental web growth machine, 

Toward the end of this report period, the new experimental web 
growth machine was placed in service and operated on a regular basis. 

The facility functioned as intended, producing web in the first run, 

Part of the melt level sensing system had not yet been installed when 
these initial nins were made, so molt replenishment was accomplished 
with manually set feed rates. Replenishment at full replacement rates 
was carried out over the period of the work day without difficulty. 

3, 2. 1.4 Programmed Start of Growth 

A controlled and reproducible start-up sequence Is desirable for 
reliabio, routine web growth. It Is also important to be able to iso- 
late this variable in studies of crystal quality as a function of other 
growth parameters. 

Typically, the growth Initiation Involves several minutes of 
coupled changes in melt temperature and growth velocity, as suggested 
schematically in Figure 20. The melt is first undercooled below the 
equilibration or "hold" temperature T^ to produce a button of suitable 
width. Tlie growth velocity is Increased to thin the button and start 
web propagation. Both temperature and velocity are then adjusted to 
their final steady state values. Small variations in the sequencing and 
rates of change in these manipulations can make the difference between 
"good" starts and "poor" starts. Thus, the quality and reproducibility 
of starts when performed manually is dependent on the skill of the 
operator. 

Programmed start of the growth operatlon^^^ has the potential 
for several important advantages compared to manual growth starts, 
including (1) repeatable wide starts, (2) reproducible high crystalline 
quality, and (3) loss stringent operator skill and training. To 
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Pull Speed Temp. Control Set Point 


Curve 725246-A 



Figure 20 — Typical Program Sequence for Web Growth Initiation 
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evaluate programmed etart of growth, we have iaterfaccd a Leedo ai^d 
Northrop 1300 Proceas Programmer with the speed and temperature control 
circuits of a wob growth furnace* This programmer Coaturos two synchro- 
nized programs which are independently set to adjust pull speed and 
temperature. 

The system has been used to provide high-quality starts of 
growth with no operator input. The interfacing has been designed to 
allow synchronized ("bumpless'') transfer from "manual" to "programmed" 
and back to "manual" status, thus permitting a single programmer to be 
transferred from furnace to furnace as required. 

The circuits designed to accomplish the "bumpless transfer" 
between manual and programmed modes are depicted in Figures 21 and 22. 

In actual operation, the transfer from the remote or programmed mode to 
local, manual control Is achieved by setting the y-jltage at V In each 
case to zero with the selector switch in remote. For temperature 
control, this is done by adjusting the alternate set-point potentiometer 
to match the output of the local set-point potentiometer. In the cose 
of pull speed, the local control potentiometer is adjusted to match its 
output to that of the interface amplifier. The locol/reraote switches 
can then be set to local and the crystal growth maintained in the 
conventional manner throughout the rest of the run. 

3. 2. 1.5 Long-Term Growth Stability 

Unless compensated in some manner, a very gradual long-term 
temperature drift may occur daring web growth. Without compensation, 
this temperature drift would eventually cause interruption of the web 
growth process, necessitating occasional restarting of the web crystal. 
Long-term drift occurs for at least two reasons; (1) gradual change of 
the temperature sensor calibration and (2) gradual change of the thermal 
properties of the system, notably the emlssivity of heat shields. 
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Figure 21 — Temperature Control Transfer Circuit 
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Figure 22 — Full Speesj Transfer Circuit 
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Wo have toe tod two raeUhods to oenao and Gompenaote the ays tom 
for oraall thermal drlfta, monitoring dendrite thlckneao ond monitoring 
wob width. Both mothoda work although we believe the latter approach 
of fora aorac diotlnct advantageo for control. 

Dendrite thickneao monitoring by vlaual oboorvatlon haa been in 
uae for quite a while aa a guide to adjuot growth and ia a oenaitive 
indicator of the thermal conditiona in the ayatera. With the unaided eye 
or optical equipment, the dendrite la viewed cdge-on; a thickening 
dendrite Irapllea the liquid ia too cold, a thin dendrite that it la too 
hot. Temperature adjuatmont con readily be made to corapenaato for drift 
due to aonaor calibration change or heat ahield einiaalvity variation 
with time. We recently attempted to quantify theae obaervations by 
employing the use of a linear array of diode aenaors coupled with an 
optical system to measure dendrite thickness. 

EG&G RETICON manufactures such diode arrays and the associated 
circuitry and so was invited to Wcstlnghouse R&D to demonstrate the 
equipment required. This equipment included the lenses needed to focus 
a magnified imago of tlie dendrite on the diode array. 

The demonstration was successful in generating a signal directly 
related to dendrite thickness. However, the light from the glowing melt 
was not adequate to obtain a well-defined image and an external light 
source had to be used. With external back lighting of the web sample In 
the furnace, an image of the dendrite having good contrast was obtained. 
With front lighting the image was discernible but the contrast was 
barely adequate. In both coses alignment of the optical system was 
quite Important. Where back lighting was used, the orientation of the 
web was critical. 

An alternative, and we feel better, method for temperature 
compensation is to monitor and control web width. Here the web is 
viewed face-on, a geometrically more favorable technique. With a 
suitably designed lid/shield configuration, our experience has shown 
chat web width can be maintained constant within a fraction of a 
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mtXXlmcter so Xong as the temperature lo eonotant. Tliuo, web width Is a 
reliable indicator for temperature adjustments. Suitable optics can bo 
used to enhance the width measurement. 

Thuo> wo have concluded that width rather than dendrite thick- 
ness sensing is preferable for drift qomponsatlon. The method can be 
applied in the form of fully outoroabic closed-loop control or, alterna- 
tively, con bo applied in the form of infrequent operator adjustment as 
needed. 


3. 2. 1.6 Test of Feed Stock Pellets from Kayex Corporotlon 

A sample of silicon in the form of shot made from semiconductor 
grade silicon provided by Kayex Corporation. The pellets are 
roughly spherical with a diameter range of about 1.6 to 2.5 mm. In the 
first of two growth runs, the Kayex material was used for replenishment 
stock only, the original crucible charge being our standard polycry- 
stalline silicon load. Growth behavior was normal during replenishment, 
and the pellets all molted readily in the crucible feed compartment. 
Because the holes in the pellet feed were far from optimum for the size 
of the pellets and the pellets were not screened to a narrow size range, 
multiple pellet drops frequently occurred. This was no prv)blera, as the 
feed rate could be adjusted to a suitable average value, and no diffi- 
culty in melting was observed. However, by late in the day, all throe 
feed holes in the pellet feeder became jammed with three pellets each 
and replenishment ceased. This problem could be readily solved with 
some adjustment In the pellet feeder hole size and some screening of the 
feedstock, particularly to remove "doubles” (two pellets stuck together). 

In the second run, the crucible charge contained 31% of the 
Kayex pellets in ordei. to test melting characteristics. 'Hie molt down 
was clean, with no trace of surface scum, indicating that the pellets 
were clean. (The pellets were used as supplied.) 

On the basis of these preliminary tests, this material appears 
to be an excellent form for melt replenishment stock in the silicon web 
process . 
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3.3 Experimental Silicon Web Sheet Growth Machine 

I I ' • I II I I I I II ■ I III! Til —71 ■ ■■ " ■ I — ■‘T— —— T ^-mn-fr - ■ i -i it rr in my i 

A new silicon web growth machine has been built and is ready for 
use in oxperiraentoj web growth. 

3.3.1 Purpose 

The purpose of this machine is to support the experimental part 
of a program to obtain a fundamental understanding of the silicon web 
growth process. A fundamental understanding is required for evaluation 
of the technical feasibility and the ultimate commercial potential of 
the process. The program to achieve this understanding combines closely 
coordinated eoroputer thermal modeling of the process and related experi- 
mental web growth. Until now. such an invoatlgntiQn was limited to 
experimental growth under transient conditions only. This new machine 
brings together, for the first time in a single machine, all of the 
growth control features developed by this program. With those features, 
evaluation of the fundamentals of the process can be performed under 
sustained steady state rather than transient conditions of growth. 

3.3.2 General Description of Machine 

The general configuration and placement of basic components of 
this machine, shown schematically in Figures 23 and 24, Is not unlike 
that of web growth machines developed earlier In this program. These 
general features are described in the following paragraphs. New 
capabilities required for sustained steady state web growth, which are 
unique to this machine, will be discussed in a subsequent section 
(3.3.3). 

As with earlier web growth apparatus developed under this 
program, web growth occurs within a water-cooled chamber (see Figures 23 
and 24) which contains a protective (argon) atmosphere. In this chamber, 
the silicon melt is contained within a shallow quartz crucible which is, 
in turn, contained within an induction-heated molybdenum susceptor. Ten 
kilohertz induction heating is used principally because it permits 
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Figure 24 — Experimental Web Growth Furnace — Front View 
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freedom In she use of thermul shields with minimal induction coupling 
and an aboonco of voltage breakdown and corona* Temperature profiles 
within the melt and within the growing web are determined by a system of 
thermal shields which surround both the susceptor and the first several 
centimeters of growing web. An XYZ mechanical positioner for the 
induction heating work coll is situated below the growth chamber in the 
lower part of the furnace frame. Tliio positioner is used to achieve 
f ront“to-back and slde-to-side thermal asymmetry within the growth 
system, and to achieve the desired degree of dendrite penetration into 
the melt. These capabilities arc necessary for experimental web growth 
with changing thermal configurations. 

As the web grows it is withdrawn from the furnace through a duct 
into the room atmosphere. The geometry of the duct is essentially of a 
large length-to-area ratio such that only a small exhaust flow of argon 
Is sufficient to prevent entry of room atmosphere into the growth chamber. 

As the growing web exits the growth chamber and duct, it passes 
through positioning guides and winds onto the web withdrawal and storage 
reel. The linear speed oC web growth is decermlned by the rotational 
speed of the reel, which is driven by a dlrect-current servo motor. 

The speed and dimensional capacity of this machine is more than 
adequate to meet the area rate-of-growth goal for this program, as had 
been determined by economic analysis of the silicon web growth process. 

A photograph of the machine nearing completion of assembly is show in 
Figure 25. 

3.3.3 Required Features for Sustained Steady State Web Growth 

The required feiitures, applied simultaneously, permit web growth 
to attain a steady state condition of growth for periods from several 
hours CO as many as severfil days. A growth period of several hours 
under steady state conditions is adequate for much of the fundamental 
study now in progress. A period of days, however, is required in order 
to appraise important long-term growth effects. This is important 
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because economic analysis has shown that a growth period of days is 
required in order for the process to become commercially competitive on 
a large scale. 

The required features for sustained steady state growth ate 
discussed in the following subsections. 

3. 3. 3.1 Constant Melt Level 

A basic requirement for sustained steady state web growth is to 
maintain constant temperature profiles within the melt and within the 
growing web. Because both profiles vary with melt height, it follows 
that the melt level must be held constant during web growth. Inasmuch 
as the silicon melt is continuously depleted as a consequence of web 
growth, it is necessary to replenish the melt continuously at a rate 
equal to the rate of depletion. 

In the new web growtli machine, melt replenishment is provided by 
a method shown schematically in Figure 26. In this arrangement, poly- 
silicon pellets ranging in size up to three millimeters are fed into a 
partitioned crucible consisting of a small compartment into which 
pellets are fed and a larger compartment from which the web is grown. 

The barrier which separates the two compartments has an opening below 
the liquid surface which permits liquid flow and allows tlie liquid level 
to equilibrate in the two compartments continuously. The barrier does, 
however, prevent movement of unmelted pellets to the growth compartment 
where interference with web growth would otherwise occur. The barrier 
also prevents ripples from reaching the web growth Interface as pellets 
are dropped into the melt replenishment compartment. Polyslllcon 
pellets such as those (Figure 27) prepared by Kayex Corporation for a 
parallel JPL program have been found to work quite well with this melt 
replenishment system. 

In order to maintain a constant melt level through use of the 
melt replenishment system, It is necessary to control the rate of pellet 
feeding such that it equals the rate of web withdrawal. To attain this 
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Figure 26 Simplified Sketch of Melt Replenishment System 
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contcol a closed-loop melt level sensing and control system (Figure 28) 
has been developed Cor the experimental web growth machine. 

For melt level sensing, this system utilizes a two milliwatt 
helium-neon laser mounted such that the light beam reflects from the 
melt onto a commercially available solid state position detector. Tlte 
control loop la closed with circuitry (Figures 28 and 29) which measures 
the detector output and provides continuously variable speed control, 
this system controls the melt level to within 0.1 millimeter, the value 
which has been determined as the accuracy needed for steady state web 
growth. 


3. 3. 3. 2 Constant Temperature 

Another requirement for maintaining constant temperature 
profiles (see Section 1.3.1) is to maintain a constant temperature at an 
appropriately selected location within the thermal system. The tempera- 
ture sensor used for control is located within the wall of the 
induction-heated susceptor. Temperature control instrumentation it- this 
application, similar to instrumentation used for low-frequency 
induction-heated Czochralski crystal growth, is shown in simplified 
schematic form in Figure 30. For this machine a fixed rather than 
adjustable range and zero offset are used. Temperature control accuracy 
obtained by this system is considerably better than the O.l^C range 
required by the web growth process for steady state growth. 

3. 3. 3. 3 Constant Width of Growth 

For steady state web growth, the web output must be maintained 
constant in width, thickness, and speed of growth. Control of web 
growth width in the experimental web growth machine is attained by way 
of appropriate thermal design used in combination with the constant melt 
level feature described in Section 1.3.1 and the constant temperature 
feature described In Section 1.3.2. The thermal design maintains width 
within one millimeter and is developed for widths up to three centimeters. 
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The design will be modified to accommodate widths greater than three 
centimeters. 

3.3. 3. 4 Constant Web Thickness and Speed of Growth 

For any web growth thermal configuration! the growth thickness is 
a function of growth speed. This relationship will remain fixed for 
stoody state conditions of web growth wherein the melt level and 
temperature are held conotani,. Tlius, if growth speed is maintained at a 
constant value the web thickness will, consequently, also remain constant. 
Tl\e dependence of growth thickness upon melt replenishment (with constant 
temperature and melt level) is shown in Figure 31. For the experimental 
web growth machine, a dc servo motor is used to provide precisely constant 
reel drive and growth speed. 

3. 3.3 *5 Programmed Start of Web Growth 

Programmed start of growth is not a mandatory requirement for 
sustained solid state web growth* Manually started growth can and 
frequently does result in a high-quality web crystal. Manual starting 
nonetheless does have serious shortcomings in that it requires compara- 
tively high skill by the web growth operator and, regardless of the 
operator's skill, frequently results in crystal imperfections which may 
ultimately lead to deterioration or failure of the web. Programmed start 
of growth has been found to eliminate essentially all of the manual start 
objections by producing good-quality, repeatable starts. The programmable 
start is a capability not yet fully utilized in that programs for highest 
quality and maximum width have not yet been developed. In the experi- 
mental web growth machine a two-channel commercial microprocessor-type 
progrmnmer has been interfaced with the temperature and growth speed 
control systems. 
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3.3*4 Suramotfy 

An expeclinonfcal oillcon web gcowth machine has boon built and is 
now capable of satisfying all of the anticipated experimental needs of the 
on-going program to attain a fundamental understanding of the web growth 
process under conditions of sustained steady state growth. Such funda- 
mental understanding is required in order that the process become eoraraor- 
clally competitive. Experimental web growth with this machine will be 
closely coordinated with a parallel effort which models the web growth 
process and provides guidance to experimental work. 

This machine is the first which combines all of the web growth 
Improvements generated by this program. Each major component of the 
machine has been individually and independently proven prior to inclusion 
in the new machine. Coordinated experimental web growth will begin in the 
next reporting period. 
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SILICON WEB ECONOMIC ANALYSIS UPDATE 


The basic input and assumptions of the IPEG economic analysis are 
unchanged. It has been deemed appropriate, however, to update the 
analysis in terms of comparatively -vilnor cost corrections for crucibles, 
labor, and electric power. The adjusted costs shown are in 1980 dollars 
and Include the following IPEG coefficients; EQPT-0.54, SQFT-110.6, DLAB- 
2.8, MATS-1.2, and UTIL-1.2. The cost breakdown for a three-day growth 
cycle is: 


EQPT » .54 X $15,600 $3424 

SQFT - 110.6 X 30 3318 

DLAB “ 2.8 X 4 hrs/run x $6. 30/hr x 117 runs/yr 8256 

MATS =*1.2 [$20+(68 hrs/run x $. 06/hr)] x 117 3381 

UTIL “ 1.2 X 3kW X 68 X $.05/kWh x 117 1432 

Annual cost per furnace $24,811 

QUAN (annual output per furnace) 1153M 

Area cost $24,811/1153 21.52/M^ 

Value added cost =* $21.52/(1000 x .15 x .95) I5.li/Wpk 

Polysilicon cost at $14.00/kG 3.9<}i/Wpk 

TOTAI, SHEET COST 19.04:/Wpk 


This projection of 19 cents per peak watt is 1.7'ji/Wpk higher than 
the previous projection but well below the JPL 1986 goal of 22.4t-/Wpk. 


5. SUMMARY AND CONCLUSIONS 


5.1 Summary 

The thermal stvooi? model haa been developed to the point where it 
can be applied to the design of experimental growth configurations. Wlien 
applied to previously tested configurations, Che model gave reoults in 
excellent agreement with experiment and is presently being utilized to 
design new low-stress lid and shield configurations which will be tested 
experimentally* 

The value of the modeling effort cannot be over emphasized. 
Experiraenta have shown that in some cases, small differences in lid and 
shield configuration can have a large Impact on both stress-induced 
deformation and residual stress in the web. Correlations between these 
effects and the results of calculations will enable us to concentrate 
experimental efforts on designs which the modeling results show to have a 
high probability for success. 

The new experimental web grovjth machine has been completed. This 
facility contains all the functions necessary for long-duration, steady 
state web growth, including melt level control. Thus, growth experiments 
can be carried out with fixed thermal conditions. This facility has 
undergone preliminary testing and performed as designed, producing web on 
the first run. On the basis of these tests, some minor modifications in 
the melt level control circuitry v^ere deemed desirable and these are being 
made. This facility, with the Inclusion of low-stress, width-controlled 
lid designs, will enable us to carry out growth experiments under steady 
state thermal conditions controlled to a degree not previously possible. 
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5 . 2 Future Work 

The iherraal model foi‘ cdlculaclng temperature profiles In the wob 
will be modified to incorporate more complex lid geometries and any number 
and shape of top shields. This will allow more complex configurations to 
be modeled accurately, how-stress lid and shield designs will be 
developed, adopted for width control, and tested for long-duration steady 
state growth. This sequence will progress to wider web crystals grown at 
higher speeds as designs are improved. 
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6 . NEW TECHNOLOGY 

1. Pcogramraed growth Initiation of silicon web crystals, JPL New 
Technology items SC1293 and SC1295* 

2. Thinner lids and Improved shield configurations for greater 
throughput, JPL New Technology item SC1267. 
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